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ABSTRACT 



Context. Thermal X-ray spectra of clusters of galaxies and other sources are commonly calculated assuming Maxwellian electron distributions. 
There are situations where this approximation is not valid, for instance near interfaces of hot and cold gas and near shocks. 
Aims. The presence of non-thermal electrons affects the X-ray spectrum. To study the role of these electrons in clusters and other environments, 
an efficient algorithm to calculate the X-ray spectra is needed. 

Methods. We approximate an arbitrary electron distribution by the sum of Maxwellian components. The decomposition is done using either 
a genetic algorithm or an analytical approximation. The X-ray spectrum is then evaluated using a linear combination of those Maxwellian 
components. 

Results. Our method is fast and leads to an accurate evaluation of the spectrum. The use of Maxwellian components allows to use the standard 
coUisional rates that are available in plasma codes such as SPEX. We give an example of a spectrum for the supra-thermal electron distribution 
behind a shock in a cluster of galaxies. The relative intensities of the satellite lines in such a spectrum are sensitive to the presence of the 
supra-thermal electrons. These lines can only be investigated with high spectral resolution. We show that the instruments on future missions 
like Astro-H and IXO will be able to demonstrate the presence or absence of these supra-thermal electrons. 

Key words. Acceleration of particles - Radiation mechanisms: general - X-rays: galaxies: clusters - X-rays: general 
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> 1. Introduction 



Most of the visible baryonic matter in clusters of galaxies is 
in the form of a hot, tenuous gas. Usually it is assumed that 
this gas is in collisional ionisation equilibrium. In the cen- 
tral parts the density is high enough and the timescales long 
enough to fulfil the conditions for collisional ionisation equi- 
librium. In the outer parts this may not always be the case. 
Freshly accreted gas may still be ionising. For a review o n equi- 
librati on processes in such tenuous plasmas see Bvkov et al. 
(120081) ■ Nevertheless, usually still a Maxwellian electron dis- 
tribution is assumed to be valid. However, there are situations 
where this is not an obvious constraint. For instance, in cool- 
ing cores of clusters, electrons from hot gas may penetrate the 
colder parts. Also, when shocks are present, deviations from 
a Maxwellian distribution may occur associated with the tem- 
perature gradients or due to particles accelerated by the shock. 
Shocks can be due to merger activity, AGN jets, and in cluster 
outskirts accretion shocks may occur. As the Coulomb ther- 
mal relaxation times increases with electron energy E as 
supra-thermal electrons are relatively long-lived and may yield 



a pressure contribution that is potentially interesting for mass 
profiles. The presence of such supra-thermal electrons can be 
revealed by excess emission in satellite lines. 

In this paper we describe a way to model the emerging X- 
ray emission spectra for the case of a non-Maxwellian elec- 
tron distribution. Although our main focus is here on clusters 
of galaxies, the same procedure can be applied to different cir- 
cumstances, for instance in stellar coronae, superno va remnants 



and th e G alactic ridge; for that last case see also iMasai et al 
(I2OO2I) and lTatischefH (l2003h . 



2. X-ray spectra for supra-thermal electron 
distributions 

To calculate an X-ray spectrum, essentially two steps are 
needed: 1) determine the ionisation balance and 2) calculate the 
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corresponding spectrum. For more details see e.g. 'Kaastra et al 
(2008.) . 

We first consider the ionisation balance. In order to deter- 
mine that, it is necessary to evaluate the collisional ionisation 
and recombination rates for all ions. These rates are derived by 
integrating the relevant cross sections over the electron distri- 
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bution. For instance, the coUisional ionisation rate Cdi is given 
by 



Cr 



OO 



(E)vf{v)du, 



(1) 



where «e and rii are the densities of electrons and the ion / con- 
sidered, cri(E) is the energy-dependent collisional ionisation 
cross-section, E = m^v^ 12 is the kinetic energy of an electron 
colliding with the ion and /(i>) is the electron velocity distribu- 
tion. 

In the most simple approximation, the cross section for col- 
lisional ionisation can be written as (lLotj|l967 ): 



CTiiE) = 



flHs in{E/I) 
£7 ' 



(2) 



where E is the kinetic energy of the free electron, Hs is the num- 
ber of electrons in a given atomic shell and the normalisation 
a - 4.5 X lO"^'* m^keV^. The cross section (|2|i can be inserted 
into ([T]), and the integration can be done analytically for sim- 
ple electron distribution functions Uke Maxwellians or power 
laws, resulting in an expression proportional to £'i(//A:r)//r"^, 
where / is the ionisation potential of the relevant shell, and 
Ei(x) - f 'e 'df is the exponential integral. Eq. (|2| is ob- 
viously too simple in realistic cases. But more sophisticated 
approximatio ns to the c ollisio nal ionisation cross section can 
be made, e.g. lYoungei ( 1981 ). that still allow analytical inte- 
gration for Maxwellian or power law electron distributions: 

m/VK^) = A(l - l/u) + B(l - l/uf + Clnu + Dlnu/u, (3) 

with u = E/I and A, B, C and D adjustable parameters. The use 
of analytical integration has the obvious advantage of enhanced 
computational speed, as during the evaluation of a spectrum 
many atomic shells of hundreds of ions need to be taken into 
account. 

A similar treatment can be made for other relevant rates, 
like the collisional excitation rates that are important for the 
line emission. However, some rates cannot be written in a form 
that allows analytical integration over a Maxwellian electron 
distribution. A good example are the radiative recombination 
rates. We recall that the cross section for radiative recombina- 
tion (Tbf can be expressed in terms of the photoionisation cross 
section through Milne's relation as 



(Thfiv) = £'icrti,(£'ph)/(TOeC^'«ef^), 



(4) 



where crfb(£'ph) is the photoionisation cross section at the pho- 
ton energy Ep^ - I + E. Although there exist analytical ap- 
proximations to th e photoionisation cross sections (for instance 
Verner & Yakovlevnl995) . they are too complicated to allow for 
analytical integration of (|4]l over a Maxwellian electron distri- 
bution, in particular as they are more complicated than a simple 
power law near the ionisation edges (Cooper minima). Even for 
the hydrogenic case, where the cross section can be calculated 
analytically, its shape is too complex for analytical convolution 
with the electron distribution. Therefore, the integration must 
be done numerically. In practice, often the exact cross sections 



(including also resonances near the edge) are convolved numer- 
ically with a Maxwell distribution, and the resulting rates are 
then approximated by analytical functions. 

It is clear that we need to follow an alternative approach 
for non-Maxwellian electron distributions. There have been 
attempts to make a generalisation of the Maxwell distribu- 
tion. For instance, Porquetet al. (2001) use a combination of 
a Maxwellian at low energies with a power-law at high en- 
ergies, but they onl y consider the io nisation balance, not the 
resulting spectrum. iProkhorov et al. I 12009) use a similar ap- 
proximation for the electron distribution, but only consider 
continuum emission and the Fexxv and Fexxvi ls-2p blends. 
Owocki & Scudder (1983.) use the so-called kappa-distribution 
(r01bertetal.il967h : 



3/2 



l_nh_\ 



r(^+i) 



1 + 



E V 
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(5) 



This distribution tends to a Maxwellian for — > oo and it 
has a power law tail proportional to £ " ' for high energies. 
However, the above approach has also clear disadvantages. 
Not every nonthermal electron distribution can be cast into 
the shape of (|5]), and not all rates can be convolved analyt- 
ically with ([5|l, in partic ular again the recombination rates. 
lOwocki & Scuddeii ( 1983 ) had to assume very simple analyti- 
cal approximations to the cross sections, like Q, in order to be 
able to evaluate analytically the plasma rates, which are given 
in terms of hypergeometric functions. Therefore we propose a 
different approach. 

3. The multi-Maxwellian approach 

Plasma codes generally may contain (ten)thousands of spec- 
tral lines, and at high spectral resolution, spectra contain many 
bins. Therefore, numerical integration for each rate over the 
electron distribution is a very time consuming task, and not of 
practical use, in particular when spectral fitting is performed 
and the model needs to be evaluated many times. 

But all relevant rates in the case of the collisionally ionised 
plasmas that we consider here, are proportional to the product 
ngWi, as for each of these rates the interaction (collision) of an 
electron with an ion is the driving process. Therefore, all rates 
are linearly proportional to the electron density Wg- If we de- 
compose the electron distribution into a linear combination of 
elementary components, the total rate for any process is simply 
the sum of the rates for these individual elementary compo- 
nents. Now, for Maxwellian electron distributions all rates are 
well known and fast to evaluate analytically or using analytical 
approximations. Thus, if we decompose an arbitrary electron 
distribution f{E) into a linear combination of Maxwellians: 



f(E)^Yj^igM(E,kTd, 



(6) 



with Aj the normalisation constant for component / and 
guiE, Ti) the Maxwell distribution for temperature T,, then the 
calculation of any rate (recombination, ionisation, excitation) 
is simple and straightforward. The problem is then reduced to 
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obtaining the normalisations A, and temperatures T, for an ar- 
bitrary electron distribution f{E). 

In principle, one might solve formally the integral equation 



RE) 



Using 



OQ 

/ 



A{kT)gM{E, kT) dkT. 



gM(E,kT). 



-E/kT 



^,1/2(^7^)3/2 

O can be cast into the shape of a Laplace transform as 



^'^'f(E) 
with 



= J B(y)e-^ydy, 





(7) 



(8) 



(9) 



(10) 



The calculation of the inverse Laplace transform is not always 
trivial. For instance, when we consider a mono-energetic beam, 
f{E) is essentially a delta-function, and the formal solution is a 
sine-wave with infinite frequency, which is not very practical. 
Fortunately, in most practical cases f{E) is a smooth function 
of energy, spanning a broad range of energies. 

There are various ways to solve (|9]l- There is software avail- 
able that can do the inverse La place transform numerically (for 
example Id'Amore et al. but in practice this is difficult, 

because these algorithms only work if the left-hand-side of (|9]l 
is an analytical function. However, in practice the electron dis- 
tributions are given in tabular form, and the algorithms used 
to interpolate such tables do not produce a formal analytical 
function (the requirement that the function is infinitely differ- 
entiable is usually violated). There are two solutions to this 
problem that we elaborate in the next subsections: direct ap- 
proximation of the electron distribution by an analytical func- 
tion, or direct decomposition of the electron distribution as a 
sum of Maxwelhans. 

3.1. Approximation of the electron distribution by an 
analytical function 

In practice, apart from the most pathological cases, electron 
distributions contain a core that is close to Maxwellian, plus a 
high energy tail. For an example see Sect. |4] After some ex- 
perimentation, we found that the electron distribution of that 
example can be approximated by the following series: 



2kTo 



C()e " -I- 



k=\ 



Ck 



with X the dimensionless momentum defined by 
X = ^JE/kTo, 



(11) 



(12) 



and where is a characteristic temperature corresponding to 
the electron distribution, and the parameters q and can be 
adjusted to give the best fit. 



Eq. ( fTTI ) is an analytical function, and the inverse Laplace 
transform in this case can be done analytically, yielding 

A{kT) = S^s{kT-kTo/ao) 



+ 



kTom + k/2) \T I 



Ck /roV3+W2 



(13) 



where Fix) is the F-function and 6(x) Dirac's delta-function. 

For our example, the normalisations co - C4 are 0.755, 
0.0609, 2.54, 13.3, and 17.58, respectively; the scale factors 
flo - 04 are 0.483, 152, 6843, 57.4 and 12.4. The fit is not per- 
fect, but better than 0.5% for x < 30, and better than 1% for 
X < 700. The relative difference to the exact distribution are 
shown in Fig.|2] 

3.2. Direct decomposition of the electron distributions 
as the sum of Maxwellians 

We seek an approximation to /(£) that can be written as a lin- 
ear combination of n Maxwellians, with n a given number. The 
solution can be represented as a set of pairs (r,,A,) with Ti 
the temperatures and A,- the normalisations. We will choose for 
convenience a set of increasing temperatures (Tj+i > Ti), and 
we will demand that A, > (physically allowed components). 
Most realistic electron distributions can be represented as a 
Maxwellian with high-energy tails. Therefore, the first temper- 
ature T\ should be close to the temperature Tq of this dominant 
Maxwellian component. 

This then leads to the following approach. We define a 
set of n logarithmically spaced temperature intervals (Ti ,, Ta,,) 
with T\ \ equal to Tq and for all /, we take r2,, = sTi^i and 
7^1./ = Ti^i-x - For the parameter i > 1 we adopt a value of 1.5. 
The total temperature range spanned by these n contiguous in- 
tervals is therefore To - s"Tq. For each interval, we choose an 
arbitrary temperature T, within this interval (Ju < Ti < T2,i) 
and an arbitrary normalisation A, > 0. We scale the A, in such 
a way that their sum corresponds to the total electron density. 
We then check how close the resulting electron distribution is 
to the electron distribution f{E) that we want to approximate. 
By varying all values T, within their allowed ranges (the above 
mentioned intervals) and also the values for A,, we try to get 
the best possible approximation. 

We use the _genetic algorithm PIKAIA developed by 



Charbonneaul (119951) to find the best solution. This algorithm 
needs a formal function of the parameters (T,, A,) that needs to 
be maximised. For this function we choose the inverse of the 
sum of the squared relative deviations of the approximation to 
the true electron distribution f{E). 

After some experimenting, we obtained satisfactory con- 
vergence with the following initial parameters: number of indi- 
viduals 500, number of generations 500. All other parameters 
in PIKAIA were kept to their default values. The number of in- 
dividuals represents here the number of sets (T,, A,) from which 
the iteration starts; the number of generations is the number of 
iterations during which the individual sets evolve towards the 
optimum solution. 
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Fig. 1. Electron distributions immediately after a shock with 
different Mach numbers, in the downstream region. The dis- 
tributions are expressed here in terms of the dimensionless 
momentum (normalised to the thermal momentum in the far 
upstream region). The distributions are normalised to integral 
unity. Solid line: M - 2.2; dashed line: M - 1.5; dash-dotted 
line: M - 1.2; dotted line: far upstream Maxwellian distribu- 
tion. 
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Fig. 3. Relative normalisations versus temperature for the 32 
Maxwellian components used in the approximation of the elec- 
tron distribution that is shown in Fig. |2l Histogram: using the 
analytical approximation ( fT3l ); Dots: using the solution solution 
from the genetic algorithm. 
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Fig. 2. Relative differences of approximations to the electron 
distribution with M - 2.2 of Fig. [T] Upper panel: analytical 
approximation ( fTTT i. see Sect. 13.11 lower panel: approximation 
using the sum of 32 Maxwellian components, see Sect. 13.21 



We have tested our algorithm on various different electron 
distributions, but give here only one illustrative example. 

4. Example of decomposition of an electron 
distribution into Maxwellians 

We have used electron distributions based on the models of 
Bvkov & Uvarov (Il999h for coUisionless MHD shocks taking 
into account particle acceleration and a shock precursor. We 
consider here a case for a pre-shock temperature of - 10^ K, 
electron density 10^ m"^, magnetic field 10 T, and shock 



size 100 kpc. In the simulation we assumed a pure Maxwellian 
distribution of electrons in the far upstream region of the flow 
and thus only the direct injection of electrons from the up- 
stream thermal pool to Fermi type shock acceleration is re- 
sponsible for the non-thermal electron population. The elec- 
tron distribution in the immediate post-shock region is shown 
in Fig.[T]for three different Mach numbers. The model distribu- 
tion was calculated taking account of both Fermi-type acceler- 
ation in a coUisionless shock and Coulomb losses of the elec- 
trons in both the upstream and downstream regions. It is nor- 
malised to the momentum po corresponding to the typical en- 
ergy fcro. Note that the supra-thermal electron distribution well 
away from the shock front differs from that shown in Fig.[T] be- 
cause the efficient Coulomb losses wash-out the non-relativistic 
supra-thermal electrons. 

There is a variant of the model where mildly relativis- 
tic electrons (with Lorentz factors > 30) comprising a pu- 
tative long-lived cosmic ray electron population in the ICM 
are re-accelerated by the MHD-shocks. The energy efficiency 
problem is alleviated in that model in comparison with the 
case of only direct particle injection from the thermal plasma. 
However, in this paper we concentrate mostly on the diagnostic 
of non-relativistic electrons and thus we consider only a direct 
injection model. 

For our example, we take the case of M = 2.2. In our de- 
composition using the genetic algorithm we have taken n = 32 
Maxwellians. The relative differences of the approximation 
compared to the exact electron distribution is shown in Fig.|2] 
The temperatures and relative normalisations of the solution 
are shown in Fig. |3] We have also used the analytical approx- 
imation ( fT3T l and binned it into similar temperature bins as the 
solutions from the genetic algorithm. 
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5. Spectrum for the supra-thermal electron 
distribution 

We have adjusted all models in the spectral fitting package 
SPEX that involve emission or absorption from a hot plasma. 
All these models now include an option to account for the pres- 
ence of supra-thermal electrons. They have an additional pa- 
rameter, which is the name of a file containing the temperatures 
and relative emission measures of the Maxwellian components. 
As SPEX does not use pre-calculated tables but calculates spec- 
tra on the fly, all relevant rates (ionisation, recombination, exci- 
tation) are simply calculated by adding the contributions from 
the Maxwellian components. Obviously, this process is done in 
two steps: first the composite multi-Maxwellian electron distri- 
bution is used to determine the ionisation balance, and using the 
resulting non-equilibrium ion concentrations, we calculate the 
X-ray spectrum for the non-equilibrium electron distribution. 
It is tacitly assumed here that we consider a time-independent, 
steady situation. For the example given in the previous section, 
only a few dozen Maxwellian components are needed. This al- 
lows fast and accurate evaluation of the spectrum, without the 
need to make simplifications to the atomic physics. 

The most important effects of non-thermal electrons on the 
spectrum is then a shift of the ionisation balance towards higher 
ionisation, the production of a non-thermal Bremsstrahlung tail 
on the continuum spectrum, and enhanced satellite line emis- 
sion. These satellite lines (for instance in the Fe-K band) can 
be detected easily with high-resolution spectrometers that will 
fly on future missions such as Astro-H and IXO. Their rele- 
vance as indicators for t he presence of non-therm al electrons 
was already indicated by lGabriel & Phillips (Il979h . 

To illustrate the effects of such a supra-thermal electron 
distribution on data, we simulated an XMM-Newton EPIC/pn 
spectrum extracted from a circular region with a radius of 1 ' 
centred on the core of a bright cluster with a 0.3 - 10 keV lu- 
minosity of Lx = 6.3 X 10^^ W within the extraction region, 
at an assumed redshift of z = 0.055. In the simulation of the 
spectrum, we assumed the above mentioned post-shock down- 
stream electron distribution for the Mach number of M = 2.2 
and pre-shock temperature of kT - 8.62 keV (10^ K). We as- 
sume a deep 100 ks observation. The resulting spectrum has 
very high statistics, which should in principle allow to detect 
any non-isothermality of the plasma. The best fit temperature 
of the simulated spectrum is 17.86 ±0.10 keV (2 x 10^ K), con- 
sistent with the post-shock temperature of the plasma given by 
the Rankine-Hugoniot jump condition. An isothermal model 
fits the data extremely well (reduced - ^ 02) and the non- 
thermal tail of the electron distribution cannot be detected in 
the spectrum. 

We also simulated a spectrum with the same input pa- 
rameters as observed during a deep 100 ks observation with 
the X-ray micro-calorimeter on the proposed International X- 
ray Observatory (IXO). We fitted the simulated spectrum with 
an isothermal model and obtained a temperature of kT - 
16.99 ± 0.03 keV, about 1 keV lower than the expected post- 
shock temperature. In Fig. H] we show the 6.9 - 7.0 keV part 
of the spectrum (rest-frame energies) which shows the Fe xxvi 
Lya lines and the Fe xxv y-satellite line. Enhanced equivalent 
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Fig. 4. Crosses: simulated 100 ks calorimeter spectrum for IXO 
(top panel) and Astro-H (bottom panel) as described in the text, 
for the supra-thermal electron distribution of Fig.[T] Solid line: 
best-fit model to a pure Maxwellian plasma, with temperature 
16.99 keV. Note the excess emission of satellite lines in the 
data, in particular the Fe xxv y'-line. 



widths of satellite lines are good indicators of non-thermal 
electrons. The satellite line in the simulated spectrum in Fig.|4] 
is clearly stronger than that predicted by the thermal model 
with a Maxwellian electron distribution. This exercise illus- 
trates, that in order to observationally reveal non-Maxwellian 
tails in the electron distributions, we will need high-resolution 
spectra obtained by future satellites with a large effective area. 

However, even before IXO, Astro-H (expected launch 
2014) will be able to detect supra-thermal electrons. We simu- 
lated the same spectrum for the same extraction region as above 
for the main instruments of Astro-H, again for 100 ks exposure 
time. The two hard X-ray telescopes detect the source up to 
~ 75 keV, and the spectrum is well approximated (x^ - 225 
for 223 degrees of freedom) by an isothermal model with mea- 
sured temperature of 17.64 ± 0.12 keV. Thus, the presence of 
such a small amount of supra-thermal electrons cannot be re- 
vealed as a hard tail, but it can be revealed in high-resolution 
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spectra. Fig.|4]shows the simulated spectrum for the Soft X-ray 
Spectrometer (SXS) of Astro-H. The excess flux at the Fe xxv 
j-satellite has a 3cr significance. Obviously, longer exposure 
times will enhance the significance. 

Interestingly, in all our simulations above, the best-fit iron 
abundance for the isothermal model is about 30 % higher than 
the actual abundance that we have put into our model spectrum 
with supra-thermal electrons. This holds also if we restrict our 
fit only to the Fe L-shell or Fe K-shell band. This is because 
the higher-energy electrons have less efficient line emissivity 
relative to the continuum, compared to lower-energy electrons. 
Without knowing the amount of supra-thermal electrons, which 
can be determined only from high-resolution spectra, this abun- 
dance bias cannot be resolved and will result into incorrect in- 
terpretations. 

6. Concluding remarks 

In practice most effort goes into finding a good decomposi- 
tion of an electron distribution into Maxwellians. We have indi- 
cated and illustrated in this paper two different methods: fits to 
simple analytical models that allow analytical inversion of the 
Laplace transform, and a direct decomposition using a genetic 
algorithm. 

Finally we note that all relevant plasma rates that are used 
in the SPEX code are calculated using non-relativistic approx- 
imations. For instance, ionisation cross sections are approxi- 
mated by analytical functions like (O that loose their validity 
for relativistic energies. Thus, for electron distributions con- 
taining a significant fraction of relativistic electrons, the results 
will be less accurate. Fortunately, in most situations this is not 
a problem. For instance, the electron distribution of Fig.[T]has 
a high-energy tail roughly proportional to p^^. Inserting this 
for instance into ([T]) and using ^ for the high-energy limit, 
shows that the integrand scales, apart from a logarithmic term, 
proportional to E^^, and therefore the highest energy electrons 
do not contribute much to the rates. Therefore, even though 
the approximations made to the decomposition of the electron 
distribution, in particular the genetic algorithm, are not always 
very accurate at high energies (see Fig.|2]), this affects the final 
spectrum to a much lesser extent. 

Only in astrophysical situations with a significan t num- 
ber of relativistic electrons our method will not apply. iBykov 
(l2002h has given an example of this in the context of super- 
nova remnants, considering only line fluorescence due to col- 
lisional ionisation. Good approximations to the relativistic col- 
lisional ionisation cross section of K-shell and L-shell elec- 



Astrophysical Observatory for and on behalf of the National 
Aeronautics and Space Administration under contract NAS8-03060. 
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